Geodesic Motions Versus Hydrodynamic Flows in a 
Gravitating Perfect Fluid: Dynamical Equivalence 

and Consequences 

K. Kleidis 1 and N. K. Spyrou 2 

Astronomy Department, Aristoteleion University of Thessaloniki, 
54006 Thessaloniki, Macedonia, Greece (Hellas) 
e-mail address: 1 kleidis@helios. astro.auth.gr , 2 spyrou@helios. astro. auth.gr 

February 7, 2008 

Abstract 

Stimulated by the methods applied for the observational determination of 
masses in the central regions of the AGNs, we examine the conditions under 
which, in the interior of a gravitating perfect fluid source, the geodesic motions 
and the general relativistic hydrodynamic flows are dynamically equivalent 
to each other. Dynamical equivalence rests on the functional similarity between 
the corresponding (covariantly expressed) differential equations of motion and is 
obtained by conformal transformations. In this case, the spaces of the solutions 
of these two kinds of motion are isomorphic. In other words, given a solution 
to the problem hydrodynamic flow in a gravitating perfect fluid, one can always 
construct a solution formally equivalent to the problem geodesic motion of a 
fluid element and vice versa. Accordingly, we show that, the observationally 
determined nuclear mass of the AGNs is being overestimated with respect to 
the real, physical one. We evaluate the corresponding mass-excess and show 
that it is not always negligible with respect to the mass of the central dark 
object, while, under circumstances, can be even larger than the rest-mass of 
the circumnuclear gas involved. 



1 Introduction 

It is generally believed that active galaxies are powered by the presence of central, 
massive black holes (Rees 1984; Blandford & Rees 1992). Many galaxies may have 
gone through an active phase and, therefore, massive black holes could be common 
in any active galactic nucleus (AGN) phenomenon (Holt et al. 1992; Antonucci 1993; 
Urry & Padovani 1995). Moreover, black holes are expected to be present in the 
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centers of many quiescent galaxies as well (Chokshi & Turner 1992). However, direct 
dynamical evidence for the presence of black holes in the central regions of individual 
galaxies is scarce. Stellar kinematical studies have provided tentative evidence for 
black holes only in a handfull of nearby galaxies, because of difficulties in spatial res- 
olution and the lack of knowledge on the exact shape of the stellar orbits (Kormendy 
& Richstone 1995; Ho 1998). Based on stellar dynamics, the mass concentrations 
in the central regions of these galaxies are determined through Doppler-shift mea- 
surements involving Keplerian motions. Accordingly, the velocity dispersion v(r) of 
various sources of radiation located at a distance r from the center, is usually assumed 
to follow the standard virial-theorem equation for circular geodesic motion (Peebles 
1972; Quinlan et al. 1995) 

v» = G^ (!) 



where M(r) is the corresponding distribution of mass inside the radius r and G is 
the Newton's gravitational constant. In fact, it was this law which led Sargent et al. 
(1978) to estimate, for the first time, the mass of the central black hole in M87, to 
be of the order 5 x 10 9 M & . 

On the other hand, recent observational data indicate that, in most of the AGNs, 
there exist gas clouds surrounding the central dark object and the associated accre- 
tion disk, on a variety of scales from a tenth of a parsec to a few hundreds parsecs (see 
e.g. Holt et al. 1992; Urry & Padovani 1995). These clouds are assumed to form a 
geometrically and optically thick torus (or warped disc), which absorbs much optical 
UV radiation and soft X-rays (Guilbert & Rees 1988; Rees 1995). Since 1993, many 
indications appeared in favour of the above idea, e.g. by HST imaging of gas and 
dust in NGC 4261 (Jaffe et al. 1993), by the mapping in nearby AGNs of the thermal 
molecular CO emission in Centaurs A (Rydbeck et al. 1993), or by the high-density 
tracer HCN in NGC 1068 (Tacconi et al. 1994) and M51 (Kohno et al. 1996). This 
circumnuclear gas seems to exist in either the molecular or the atomic phase, the lat- 
ter being due to hard X-ray sources in the centers of most radio-loud AGNs (Krolik & 
Begelman 1988; Maloney et al. 1994). In fact, to date the most dramatic evidence for 
the existence of a supermassive black hole comes from the VLBI imaging of molecular 
(H 2 0) masers in the active galaxy NGC 4258 (Greenhill et al. 1995; Miyoshi et al. 
1995). This imaging, produced by Doppler-shift measurements assuming Keplerian 
motion of the masering source, has allowed a quite accurate estimate of the nuclear 
central mass, which has been found to be a 3.6 x 10 7 M supermassive dark object, 
within 0.13 parsecs. However, this evidence (as well as many of the previ- 
ous ones) comes primarily from the motions of gas streams rather than 
stars. The gas is also subject to non-gravitational forces and may not follow ballis- 
tic trajectories. Indeed, it has been recently announced (Greenhill & Gwinn 1997) 
that the rotation curve of the active galaxy NGC 1068 (in which, maser emission has 
also been detected) is sub-Keplerian, suggesting that a non-gravitational component 
is also present in the maser emission. Therefore, as regards the measurement data 
relevant to the cental dark objects in the AGNs, there is some inherent ambiguity 
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(Rees 1995). 

In fact, the determination of the central nuclear mass is based on the Doppler- 
shifted radiation received, not from something like a test-particle star, but from the 
extented masering fluid sources. These sources move in the gravitational field of a 
(more or less) continuous gravitating source, such as the circumnuclear gas and the 
accretion disk, and from a dynamical point of view, the assumption that the motion of 
the masering source is Keplerian (or, in general, geodesic), seems to be inaccurate. On 
the other hand, a fluid volume element can be considered as more realistic constituent 
of a fluid source and hence, its use in a continuous medium should be physically 
preferable than that of a theoretical, ideal test-particle (Spyrou 1997a; 1997b; 1997c). 
So, the question is raised as to whether the hydrodynamical flow of a fluid volume 
element in a continuous gravitating source can be approximated or not to the motion 
of a test-particle (namely the geodesic motion) in the same source. 

Motivated by the above considerations, in the present article we examine the 
conditions under which the hydrodynamic flows in the interior of a bounded self- 
gravitating perfect fluid source may become dynamically equivalent to general- 
relativistic geodesic motions in this source. It is clear that, generally, in the interior 
of a continuous perfect fluid source the hydrodynamic flow differs from the geodesic 
motion. In this case, dynamical equivalence rests in the fact of the functional similar- 
ity between the corresponding (covariantly expressed) differential equations of motion 
(Spyrou 1999). Accordingly, the spaces of the solutions of these equations are iso- 
morphic. In other words, given a solution to a problem of hydrodynamic flow in a 
gravitating perfect fluid, one can always construct a solution formally equivalent to the 
problem geodesic motion of a fluid element and vice versa. The problem attacked is 
kept as general as possible, being solved within the context of the exact General Rel- 
ativity Theory (GRT). Therefore, the corresponding results could be useful in many 
astrophysical situations and not just in the context of relativistic galactic dynamics 
and/or cosmology. 



2 Geodesic Motions Versus Hydrodynamic Flows 

The ballistic motion of a test-particle in the gravitational field of a bounded, self- 
gravitating perfect fluid source, takes place along geodesic trajectories, on which the 
tangent vector of the four-velocity u % = ^ (uiU 1 = 1) obeys the standard geodesic 
equation in a coordinate frame x l (see e.g. Papapetrou 1974) 

^- = o & ^! + r >V = o (2) 

ds ds ki v ' 

where the Latin indices, denoting spacetime coordinates, admit the values 0, 1, 2, 3 and 
the Greek indices, denoting spatial coordinates, admit the values 1,2,3. In Eq. (2), 
denotes the covariant acceleration vector, s is the afflne parameter and Y l kl are the 
Christoffel symbols. On the other hand, the equations of motion of the hydrodynamic 
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flow in the interior of the same source, are written in form 



Vi = ( 3 ) 

where a semi-colon denotes covariant differentiation and T u is the stress-energy tensor 
of the perfect fluid, which is written in the standard form 

T ij = (£ + p)u l u j - pg ij (4) 

In Eq. (4), p is the isotropic pressure and £ is the overall mass-energy density of the 
perfect fluid. Both quantities are considered to transform like scalars under coordinate 
mappings (Anderson 1967). In the absence of shear and viscocity, £ is decomposed 
as follows (Fock 1959; Chandrasekhar 1965) 

£ = pc 2 + pli (5) 

where, pc 2 ^ is the rest-mass energy density of the total number of baryons included 
in the unit volume element, and pll is the corresponding specific internal energy 
density, the variations of which are related to the compressions or the expansions of 
the fluid. Combination of Eqs. (3) and (4) yields 

(£+p), j u l u j + (£+p) (u^ u j + u l u j . J )-p t j g ij = (6) 

Now, in order to study a possible dynamical equivalence of the hydrodynamic flow 
of a finite-volume element in the interior of a bounded, self-gravitating perfect fluid 
to the corresponding geodesic motions, first of all we need to determine the evolution 
of the fluid's covariant acceleration. In this case, it is convenient to study the flow 
motion on a hypersurface normal to the direction of u l . We do so, by contracting Eq. 
(6) with the projection operator 

h ik = 9ik - UiU k (7) 

for which, we obviously have 

h ik u l = (8) 

Then, provided that £ + p ^ 0, i.e. the total enthalpy of a fluid element is non- 
vanishing, the combination of Eqs. (6) and (7) yields 



Uk ' jUl = {£ + p) ~ UkU ^ 



(9) 



Eq. (9) is just one of the two equations of fluid motion, in which Eq. (3) splits in 
the case of a perfect fluid (Taub 1959). This equation is completely analogous to the 
conservation- of-momentum equation of Newtonian hydrodynamics (Ryan & Shepley 
1975). In a Riemannian space (gij- s = 0), Eq. (9) is written in the form 

Quu]^ = -J^(5 j k -u k v?) (10) 
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from which, on contracting with g mk , we obtain 

du 1 



, +T\ s u j u s 
ds 3 



= 9 



mk P,i 



(S 3 k - U k U 3 ) 



thus, resulting to 



du n 
ds 



+ TJ s u 3 u s 
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mk P, 3 



or, finally, 



du m 

ds 3 



u k u~ 



h mj P,j 



(12) 



From Eq. (12), it becomes evident that, the two kinds of motion are completely (not 
just dynamically) equivalent, provided that the pressure gradient in the three-surface 
normal to u l [which appears in rhs of Eq. (12)] is zero. This quantity measures the 
response of a particle to non-gravitational fields and, therefore, it is responsible for all 
deviations of flow lines from geodesies (see also Misner et al. 1973; Ryan & Shepley 
1975). 



3 Dynamical Equivalence Conditions 

By definition, dynamical equivalence between geodesic motions and hydrodynamic 
flows rests in the functional similarity between Eqs. (2) and (12). In this case, the 
spaces of the solutions of these equations are isomorphic. In other words, given a so- 
lution to a problem of hydrodynamic flow in a perfect fluid, one can always construct 
a solution formally equivalent to the problem geodesic motion of a fluid element and 
vice versa. If we manage to establish dynamical equivalence between Eqs. (2) and 
(12), then we will be able to extrapolate every result obtained on the basis of geodesic 
motions, to the more realistic context of general-relativistic fluid hydrodynamics. A 
direct application may concern observational data from many relativistic astrophys- 
ical (or cosmological) systems, such as the central cores of the AGNs (Kormendy & 
Richstone 1995; Rees 1995) or the existence and the nature of the dark matter (Saglia 
1996) etc. Dynamical equivalence between Eqs. (2) and (12) may be performed in 
two ways: 

1. By eliminating the rhs of Eq. (12). A direct way to do so is to assume isobaric 
motions, i.e. p = constant. However, in general, this is not always a physically 
necessary condition. 

2. By transfering the problem to a virtual self-gravitating perfect fluid, the metric 
of which is conformal to g k i- In terms of this fluid, Eq. (12) will be subsequently 
written in the form of Eq. (2). 
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The latter, is a very delicate method of obtaining dynamically equivalent descrip- 
tions in any gravitational field theory and not just GRT (e.g. see Spyrou 1976; Whitt 
1984; Cotsakis 1993). A conformal transformation (Penrose 1964; Hawking & Ellis 
1973) shrinks or stretches the entire manifold, introducing a change in the scale of all 
lengths and time, which can be different for the various world points, but is the same 
for all spatial directions and time at a given point. A conformal transformation of a 
metric is usually described by 

g kl = n 2 (x l ) g kl (13) 

for some continuous, non- vanishing, finite, real function Q(x l ). Accordingly, in what 
follows, we assume that there exists a virtual gravitating perfect fluid, producing a 
new metric tensor (gki), in terms of which we may write Eq. (12) in the form of Eq. 
(2), namely 

^ + T™u k u l = 0. (14) 

Next, we shall try to determine the relation between these two fluids. To begin 
with, we note that when Eq. (13) is applied, the following transformations hold 
(Hawking & Ellis 1973) 

rs = iM + ^W + ^-M'-n.) 

du m 1 d ( 1 \ 1 du m 

If = ndx^Kn) umu ' n + w ~di (15) 

With the aid of Eqs. (15), the geodesic equation (14) may be decomposed in terms 
of the original metric tensor (gki), as follows 

^— + r%u k u l = h mn -^{InQ) (16) 

ds kl dx nK 1 v 1 

Now, in view of the above definitions, dynamical equivalence between the two kinds 
of motion under consideration, implies 

h ran JL(^Q) = h ™n (17) 

dx nK ' S+p y ' K J 

To solve Eq. (17) with respect to Q is a very difficult task, especially when no 
further conditions are imposed. Nevertheless, there exists a quite general physical 
condition which may help us to reach at some (at least) particular solutions. It is 
the equilibrium hydrodynamics hypothesis, i.e. the constancy of entropy (S) along the 
original fluid's flow lines (Taub 1959; Misner et al. 1973) 

S m u m = (18) 

In fact, in relativistic astrophysics and cosmology we usually assume S to be a group 
invariant (constant in space) and therefore S = const. (Ryan & Shepley 1975). Hence, 
in what follows we consider that the flow motion is isentropic. 
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Adiabaticity of the hydrodynamical motions could be physically necessary, be- 
cause then both the thermodynamic and the matter content (the number of baryons) 
of a finite-volume element of perfect fluid remain constant during its motion (Chan- 
drasekhar 1983). Then, the general- relativistic analogue of the first law of thermody- 
namics is valid with respect to the original fluid, for dS = 0. In the case of a perfect 
fluid source, this law may be expressed in the form 



or, more conveniently, 



P, 



pc 



'£+P 



Inserting Eq. (20) into Eq. (17), we obtain 

d 



h r 



dx n 



In (Q) — In 



In Eq. (21), we perform the substitution 



Q{x k 



£+p 

pc 2 



= 



(19) 



(20) 



(21) 



(22) 



Then, the new metric (g mn ), in terms of the original one (g mn ), is written in the form 

2 



gmn 



E+P 



2$(x fe ) 



(23) 



and the dimensionless, real scalar function $(x fe ) satisfies the compatibility condition 



r$„ = o 



(24) 



Eq. (24) has a clear physical meaning: The projection of the vector $ ) n on a three- 
dimensional hypersurface normal to u n is zero. By virtue of Eq. (8), this condition 
implies that $ n is always collinear to u n . Accordingly, the general solution to Eq. 
(24) may be obtained in terms of the differential equation 



$ n = C(* m K 



(25) 



(in connection see Anderson 1967; Ryan & Shepley 1975) where C(^ m ) is & n arbitrary 
scalar function with dimensions of inverse length. The general solution to Eq. (25), 
may be expressed in the form 



${x m ) = J ({X m ) U n dx n + i{x T ) 



(26) 
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where no summation is assumed over n, and £(x r ) is a dimensionless scalar function 
of the coordinates x r (r ^ n). Obviously, the exact form of the general solution to 
Eq. (24) remains undetermined. This is not an unexpected result, since Eq. (24) is a 
system of four partial differential equations involving five unknown quantities, namely, 
the four components of u % (x n ) and <&(x n ). In this case, supplementary conditions 
should be imposed, for someone to deal with a well-defined problem. 

Nevertheless, there exists one particular solution to Eq. (24) which is of great 
theoretical interest, namely 

$(x n ) = C = const. (27) 

In this case, the combination of Eqs. (23) and (27) relates the physical metric tensor 
g mn , due to the gravitational field in the interior of a real perfect fluid source (in 
which the motions follow the laws of relativistic hydrodynamics), to the correponding 
metric g mn of a virtual fluid (in which the original hydrodynamic motions become 
geodesic). According to Eq. (23), the functional form of the metric attributed to the 
new gravitational source includes the overall enthalpy content (S+p) within a specific 
volume (^) of the original one (specific enthalpy). 

The family of conformal transformations (13) [or (23)] represents an algebraic 
group over the space of the Riemannian metrics. In this case, it is evident that the 
corresponding operation among the various group elements is the standard multiplica- 
tion between the real, finite, non-zero, scalar functions Q(x n ). There exists an identity 
element (Q = 1), as well as an inverse one (fi _1 ), such that Q^ 1 ■ Q = 1. As regards 
the identity element, it represents the degenerate case, in which the associated metrics 
(the real and the virtual) are identical, g mn = g mn . In our problem, this corresponds 
to an isobaric flow of the original fluid [e.g. see Eq. (12)]. For p = p = const., the 
first law of thermodynamics [Eq. (19)] is reduced to 

n + — = const. (28) 
P 

and the identity transformation, through Eqs. (23) and (27), implies 



C = - In 



c 2 p 



const. (29) 



Therefore, if the adiabatic flow of a finite- volume element in the interior of a bounded, 
self-gravitating perfect fluid is (in addition) isobaric, then the corresponding metric 
tensor is conformally invariant under the group of transformations (23), and the 
hydrodynamic flows are completely equivalent to geodesic motions. 

In what follows, we shall assume the validity of Eq. (27) and, unless otherwise 
stated, we shall ignore the constant factor e 2C arising in the formula for Q from the 
combination of Eqs. (23) and (27). This normalization may be performed through 
a redefinition of the proper length in the original metric, as s — > A = s e~ c . In 
concluding, the appropriate conformal transformation which guarantees the dynam- 
ical equivalence between hydrodynamic flows and geodesic motions in a gravitating 
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perfect fluid, is 



9mn 




9mn 



(30) 



provided that the following major assumptions are valid: 

1. We have considered a gravitating perfect fluid, i.e. isotropic in pressure, with 
no shear and viscocity. 

2. In this source, the hydrodynamic flow motions are adiabatic, i.e. S is constant 
along the fluid flow lines. 

3. The gravitational field in the interior of this fluid corresponds to a Riemannian 
spacetime, i.e. gki-m — and det\gki\ ^ 0. 

4. The corresponding metric tensor is a solution to the Einstein field equations, 
through T. fc / = 0. 

4 Reduction to the Newtonian Limit 

In many cases of particular astrophysical interest, the corresponding weak field-limit 
results seem to be more appropriate. For example, it has already been shown (Spyrou 
& Varvoglis 1982) that, in the context of relativistic galactic dynamics and astro- 
physics, the use of the first post-Newtonian approximation is more than enough in 
most of the cases. In fact, for astrophysically interesting sources with non- negligible 
relativistic characteristics (such as e.g. the nuclei of the giant elliptical galaxies) the 
small post-Newtonian expansion parameter (for the solar system being 10~ 6 ) is 
still very small (~ 10~ 3 <C 1). Accordingly, we will examine how the above results, 
obtained within the context of the exact GRT, may affect on the corresponding New- 
tonian limit. In this limit, as regards the metric tensor of the original gravitating 
perfect fluid, we have (Schiff 1960; Chandrasekhar 1965) 



where the symbol O n denotes terms of order n in the small parameter - <C 1, with 



with v denoting the measure of the spatial three-velocity vector and U the gravita- 
tional potential, related to the mass density p through the standard Poisson equation 



g 00 = I-2- + O4 



c 



90a = C 3 

9 a/3 = —$a(3 + 2 



(31) 



e 2 = max (v 2 , U, IT, -) 



(32) 



\7 2 U = -AnG p 



(33) 
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According to Eq. (30), the components of the metric tensor attributed to the 
virtual perfect fluid, may be written in terms of the corresponding components of the 
original one (31), as follows 



fi'oo = 



'S+p 



1 + 2- 



U 



1 - 2^ + 4 



n + ^ 



1 - 2^ + 4 



= 1-2- + 



(34) 



In this case, from Eq. (34) we observe that, g 00 is formally equivalent to g 00 , provided 
that U is replaced by the virtual gravitational potential 



u = u- n + 



v 



(35) 



Eq. (35) is identical to the corresponding Newtonian result of Spyrou (1997a; 1997b; 
1997c). In the same fashion, we obtain 



9oa — 3 — g 0o 



(36) 



and 



g a /3 = (Sa/3 + 2 ) 

= Sa/3 + 2 
= 9a(3 



i+ 4( n +9+° 4 



(37) 



The set of Eqs. (34) - (37) represents the appropriate weak field-limit transformation, 
relating the conformal metrics of the real and the virtual perfect fluid, under the action 
of the group (30). In analogy to Eq. (33), we may define the virtual rest-mass density 
(p), which produces the corresponding gravitational potential U, through 

V 2 U = -4nGp (38) 



which, in view of Eq. (35), yields 



5 The Conformal Stress-Energy Tensor 

The conformal transformation (30) affects both the geometry and the physical quan- 
tities relevant to the matter content of the gravitating source (the stress-energy tensor 
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of the perfect fluid), in a way that the dynamical description of gravity in the virtual 
case (the functional form of the Einstein field equations) remains invariant. In this 
respect, we consider that g mn is a solution of the conformal Einstein field equations 

Kmi - ^gmiK = -K% nl (40) 

where also 

8irG 

ft — 7~ 

C 4 

By virtue of Eqs. (40), it is evident that the corresponding conformal stress-energy 
tensor {T m i) satisfies the conservation law T~Vf = 0, where now the covariant derivative 
is taken with respect to the conformal metric, g m \. In this case, we may deduce the 
exact form of the stress-energy tensor, responsible for the gravitational field of the 
virtual fluid, in terms of the corresponding tensor of the original source. 

In four dimensions, the influence of a conformal transformation of the form (13) 
on the dynamical quantities describing the gravitational field (the Riemann tensor 
and its contractions) is (Hawking & Ellis 1973) 

^^"'"Q^+^^M (41) 

and 

TUi^ + e^nn (42) 

where the d' Alembert operator (□), with respect to the original metric, is given by 

= ( V^g 9 mn n 

With the aid of transformations (41) and (42), Eq. (40) reduces to 

r t 2 an 2 1 o 

1ml ^ml 9ml ^ 7-7 ! nd 

K \l K \l 

+ --Ln m ^ l --g ml ^g ns n n n s (43) 
and therefore, the trace of T m i is written in the form 



f = ~g ml f ml 



1 

fi2 



3c 4 an 

T H 

AtcG tt 



(44) 



Clearly, Eqs. (43) and (44) hold for an arbitrary form of both T m i and T m i- However, 
in connection to what previously stated, we may assume that the conformal stress- 
energy tensor can be written in the form of a perfect fluid, namely, 

f ml = {S + p) u m ui - p g m i (45) 
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for which we have 



f = £-3p 



(46) 



where, £ and p are the corresponding energy density and pressure. In this case, Eq. 
(45) results to 



3c 4 an 



(47) 



By virtue of Eq. (47) we perform the ansatz that, under the action of the conformal 
group (30), both £ and p transform as follows 



p 



1 

1 

f} 2 



£ + a 



c 4 an 



v- 



AnG D, _ 
3 -a, c 4 an 



Aug n 



(48) 
(49) 



where a is a numerical constant, the exact value of which can be determined in terms 
of the corresponding weak field-limit results [Eq. (39)] (see also Spyrou 1999). In 
analogy to Eq. (5), we assume that £ may be decomposed as follows 



£ = pc 2 (i + ^) 



(50) 



where pIT corresponds to the specific internal energy density, related with the con- 
tractions or the expansions of the virtual fluid. Then, in the weak field-limit, Eq. 
(48) is reduced to 



p(i + -j) = P 



i-4( n + 2 - 

& p 



+ 



+ a 



AnG 



<r p 



l_ & 

c 2 dt 2 



which, in the stationary case, to lowest order in \ yields 



p = p — a 



1 Wn + * 



AixG 



Now, by comparison of Eqs. (39) and (52), we obtain a 
(49) are finally written in the form 



n 2 



£ 



p = n 2 [ p - 



c 4 nn\ 

a^gIT) 

c 4 nn\ 
3ttG IT) 



u + p 

p. 



(51) 



(52) 

-1. Hence, Eqs. (48) and 

(53) 
(54) 
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obviously satisfying Eqs. (47) - (49). As regards a self-gravitating perfect fluid, 
Eqs. (53) and (54) represent the transformation law of the corresponding energy 
density and pressure, under the action of the conformal group (30). It is worth 
noting that, according to Eq. (54), the pressure of the virtual perfect fluid may not 
be constant (or vanish), although the motions in its interior are geodesies. This 
result is a distinguishing feature of the present analysis, with respect to what we 
have encountered so far in relativistic astrophysics and cosmology, namely, that the 
isobaric flows are geodesies (see e.g. Weinberg 1972; Misner et al. 1973; Papapetrou 
1974; Ryan & Shepley 1975; Narlikar 1983) and so, it justifies the use of the term 
"virtual", in denoting the matter content responsible for the conformal metric. 

In relation to the measurement data relevant to what we observe and determine 
observationally in the central regions of the AGNs, the theoretical result (53) has a 
clear physical interpretation: £ is what we actually measure by assuming geodesic 
motions, while £ is the real physical quantity associated to the fluid matter content 
existing in those regions. In this respect, the quantity 

c 2 an 

p * = ~4^ 7T (55) 

may be identified as an extra inertial energy density, associated to the contributions 
from the internal physical characteristics (internal motions, pressure, thermodynamic 
content) of the original fluid, to the measured quantities. In fact, pi describes the way 
the various physical characteristics of the source (beyond its rest-mass density) act 
as additional gravitational sources, thus affecting the geodesic motions (Spyrou & 
Dionysiou 1973; Spyrou 1997a). In this case, Eq. (53) is written in the form 

£ = ^ {£ + Pi<?) (56) 

and the question arises as to whether the observationally measured energy £ is being 
overestimated or not, with respect to the physical quantity £. 

According to Eqs. (55) and (56), the answer to the above question depends on the 
sign of Ofl, which can be determined by taking into account the so called weak energy 
condition. It is known (Hawking & Penrose 1970) that, for any physical system (as it 
is the case for the original perfect fluid, described by T m i), the weak energy condition 
is valid, i.e. the locally measured energy is non-negative 

T ml u m u l =£>0 (57) 

Admiting that the virtual fluid also represents a physical system, we demand that a 
condition similar to Eq. (57) is valid with respect to T m i as well, namely, 

f ml u m u l = £ > 0. (58) 

With the aid of Eq. (53), Eq. (58) is decomposed in terms of £, as follows 

c 4 an 

£ *Sg1T < 59 > 
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Now, by virtue of Eq. (57), the condition (59) is valid for all E > 0, provided that 

an < o (60) 

Therefore, we conclude that p« > and according to Eq. (56), the observationally 
determined energy (S) is being overestimated as compared to the real, physical 
one {£). The equality in Eq. (60) corresponds to the identity transformation case, 
representing an isobaric flow of the original fluid, for which, from Eqs. (20) and (22), 
we readily obtain Ofl = 0. 

6 Some Astrophysical Applications 

Recently, it has been proposed that the centers of the AGNs contain large amounts of 
gas, in the form of warped disks or tori, which absorbs much of the non-thermal power 
generated in the central core (probably due to accretion phenomena) and reradiate 
it (Antonucci & Miller 1985; Guilbert & Rees 1988; Holt et al. 1992; Jaffe et al. 
1993; Rydberg et al. 1993; Tacconi et al. 1994; Kohno et al. 1996; Greenhill & 
Gwinn 1997). Direct mapping of these gaseous structures can be possible with high 
spatial resolution spectroscopic observations of water masers (Urry & Padovani 1995). 
Indeed, VLBI imaging of water masers in NGC 4258 has provided a firm evidence 
for rotating molecular gas in a region less than 0.25 pc from the center of the galaxy 
(Greenhill et al. 1995; Miyoshi et al. 1995). Assuming Keplerian motion of the 
masering source, this imaging has allowed a quite accurate estimate of the central 
nuclear mass, which has been found to be a 3.6 x 10 7 M supermassive dark object, 
within 0.13 pc. However, combining these observations with the theoretical results 
obtained in the present article, it becomes evident that assuming Keplerian (geodesic) 
motions in the central region of the AGNs corresponds to work within the context 
of the virtual fluid. In this case, the virtual results can be extrapolated to the real, 
physical ones (based on hydrodynamic motions) through Eqs. (30), (53) and (54). 

As a direct application, we will calculate the mass-excess rrii, arising in the de- 
termination of the central nuclear mass of the active galaxy NGC 4258, due to the 
assumption of Keplerian motion of the masering source. In fact, represents the 
contribution of the fluid's internal characteristics (p^ to the measured quantities. 
The corresponding calculations will be carried out in the weak field limit, since the 
various measurements based on observational data relevant to the central region of 
the AGNs, are usually derived by methods of Newtonian gravity. In this limit, the 
superposition principle is valid and therefore, each gravitating source (the black hole 
and/or the circumnuclear gas) is expected to contribute separately. For the sake of 
simplicity we shall assume that the nuclear region is spherically symmetric, so that 
all the physical quantities describing the galaxy's interior are functions of the radial 
coordinate r. We consider an idealized model, according to which the circumnuclear 
gas surrounding the central dark object can be represented by a stationary perfect 
fluid which undergoes adiabatic, hydrodynamical flow under the influence of its own 
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gravitational field, plus the gravitational field of the black hole (Novikov & Thorne 
1973). In this case, dynamical equivalence between hydrodynamic flow and the Kep- 
lerian motion of a fluid's volume element implies that Eq. (35) holds, provided that 
the original gravitational potential U is now replaced by Ubh + Uq, where 

u BH = G Mm. (6i) 

r 

is the gravitational potential due to a (Schwarzschild) black hole of mass Mbh, and 
Uq is the corresponding potential of the self-gravitating gaseous matter in this region. 
Now, using Poisson's equation, Eq. (39) takes on the form 

' p = p+ ^G V2{Ii + l p ) -h MB « V2{ l ) (62) 

where V 2 (±) = -47r<J(r). The integration of Eq. (62) over a spatial volume V of 
linear dimension 0.25 pc, will give us the virtual rest-mass (m) of the central region 
in NGC 4258, namely, 

m = (m + M BH ) + TOj (63) 

where m is the total rest-mass of the gas included in this region, while m; represents 
the relative error introduced in the determination of the central mass by ignoring 
the contribution of the fluid's internal characteristics. Therefore, it becomes evident 
that the measured mass m is being overestimated with respect to the real, physical 
quantity m + Mbh- 

To determine the integral rrii = J v pid 3 x, we consider that the lower limit of 
integration corresponds to the radius of the innermost stable circular geodesic 

6GM BH (aA , 
r = 3R S = = (64) 

cr 

We see that the value of r is directly proportional to the mass of the central black 
hole. In general, the estimated masses of the galactic nuclear dark objects fall in the 
range 10 7 — 10 9 M & (Kormendy & Richstone 1995). Accordingly, the corresponding 
values of r range from 10" 5 pc to 1CT 3 pc (in connection see Holt et al. 1992). In 
the present article we adopt a mean value of 1CT 4 pc. As regards to the upper limit 
of integration, we take r max = 0.25 pc, which corresponds to the outer boundary of 
the masering annulus in NGC 4258 (Miyoshi et al. 1995). Moreover, by virtue of the 
adiabaticity condition (19), we obtain 

' V 2 (U + P)= Pl = -^-V-(-Vp) (65) 



4nG y p' r AixG p 

In the special case of adiabatic motions, the pressure is related to the rest-mass density 
through the equation of state 

p = k p 1 (66) 
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where k > and 7 > 1 are constants. Now Eq. (65) takes on the form 



P * = ^G P 



V 2 p+(7-2)-(Vp-Vp) 



(67) 



In particular, we consider that the proper-mass density distribution can be described 
by a Plummer-type function (Maoz 1995) of the form 

p = p (l+x 2 )-% (68) 

where x — ^, while p and n are possitive constants. In this case, Eq. (67) results to 

Aryn p^p 1 1 



and the associated mass is written in the form 



{[n( 7 -l) + l]x 2 - l} 



rrii 



k^n 

~G~ r ° Po 



1 p 

Jxq 



x 2 {[n( 7 -l) + l]g; 2 - 1} 

(l + a ;2)S(7-l)+2 



dx 



(69) 



(70) 



where x — 1 and x max = 0.25/10 4 = 2500. The integral on the rhs of Eq. (70) can 
be expressed in terms of elementary functions only if 



n 



(7 — 1) + 2 = £ : an integer 



(71) 



(Gradshteyn & Ryzhik 1965) and in fact, on physical grounds (n > and 7 > 1), we 
must have £ > 2. The case £ = 2 corresponds, for every value of n, to an isothermal 
flow of the perfect fluid, i.e. 7 = 1. We shall consider this case separately. For £ > 2, 
the general solution to Eq. (70) is written in the form (Gradshteyn & Ryzhik 1965) 



k'yn ! 
77- r ° Po 



2£-3 



4(£-l) 



x 



2£-5 (1 + x 



2\£-l 



(2£-3)(2£-5) (l + x 2 y- 1 



+ 



+ 



4(€-l) 



g(2£-l)...(2£-2j + l) 1 



(2€-l)(2€-3)(2€-5) ^ {t-l)...{l- 3) 
4(£-l) (2£-3)!! 1 



2j (1 + x 2 y-j 



(2£-3)(2£-5) (^-1)! 2*- 1 



tan 1 x 



(72) 



.'-•0 



Clearly, Eq. (72) is too complicated to be useful in astrophysical applications. On the 
other hand, the isothermal case (£ = 2) seems to be more appropriate. Indeed, recent 
results regarding the properties of mass accretion in NGC 4258 (Neufeld & Maloney 
1995), indicate that the adiabatic flow of the the gaseous fluid is actually isothermal, 
with a constant speed of sound c s = 7 km/ sec. In this case, 7 = 1 and, therefore, Eq. 
(72) is written as 



rrii 



kn 

^r 



,x 2 + 2 



x 



x 2 + 



-) — 2tan 1 x 



(73) 



X'() 



16 



which, by virtue of Eq. (64), results to 



rm = (jfen).1.66 x 10~ 17 M BH (74) 

However, in the case of isothermal flow, we furthermore have 

k = c 2 s (75) 

and, therefore, as regards NGC 4258, adopting the typical values c s = 7 km/ sec 
(Neufeld & Maloney 1995) and n = 5 (Miyoshi et al. 1995) we finally obtain 

m 4 ~5x 1(T 5 M BH (76) 

This very small relative error, is consistent with observations and in fact, provides a 
theoretical explanation for the almost perfect Keplerian rotation-curve observed for 
the gas, in the central region of this particular galaxy (Greenhill et al. 1995; Maoz 
1995; Miyoshi et al. 1995). Nevertheless, had we applied the same model in NGC 
1068, where the outer radius of the masering source is r max = 1 pc (Greenhill & 
Gwinn 1997), we would have obtained 

m 4 ~2x 1(T 4 M BH (77) 

while, in the extreme case of NGC 4261, where (non-masering) circumnuclear gas and 
dust appears to extend up to 

r max ~ 150 pc from the center (Jaffe et al. 1993), we 
would have accordingly obtained 

m } ~3x 10~ 2 M BH (78) 

Therefore, is not always negligible compared to the mass of the central 
dark object and it can account from a few hundredths of thousanths to several 
hundredths of Mbh, depending on the linear dimensions of the circumnuclear gas 
observed in the AGNs. 

On the other hand, the assumption of Keplerian motions in the cental region of 
the AGNs is compatible with the existence of circumnuclear gas around a black hole, 
only if the central dark object contains at least 98% of the galactic nuclear mass 
(Greenhill et al. 1995). In fact, calculations based on Newtonian dynamics indicate 
that the total rest-mass of the gas is m ~ ^ Mbh (Maoz 1995). In this case, we may 
also calculate the ratio ^ for the particular galaxies considered above. Accordingly, 
for NGC 4258, we find 

mi = 2.5 x 10~ 3 m (79) 

while for NGC 1068 and NGC 4261 we obtain mj ~ 10~ 2 m and ~ 1.5 m, re- 
spectively. Therefore, m« can be comparable to or even larger than the total 
rest-mass of the circumnuclear gas involved. 

At the outcome, we have to point out that the above idealized scheme is rather 
naive in realistic situations. In fact, it is expected that in the central regions of the 



17 



AGNs the flow motion is not adiabatic at all. Significant energy losses may occur, due 
to thermal bremsstrahlung, synchrotron radiation and/or radiative transfer (Novikov 

6 Thorne 1973), leading to a loss of angular momentum, so that, ultimately, the 
matter surrounding the central dark object will end up into the hole (Shakura & 
Sunyaev 1973). It has been found (Neufeld & Maloney 1995) that, for a mass accretion 
rate of the order rh = 7 x 10 -5 M yr^ 1 , material located at a distance 0.25 pc from 
the center will fall into the hole after 4 x 10 6 yrs. Therefore, the case of a viscous 
accretion disk seems to be more appropriate than the spherically symmetric perfect 
fluid distribution, but it will not be considered here. 

7 Discussion and Conclusions 

Although recent developments in galactic dynamics indicate that the center of the 
AGNs contains large amounts of gas (Urri & Padovani 1995), the observational de- 
termination of the central nuclear masses is based on the assumption of pure geodesic 
motions (Kormendy & Richstone 1995). Stimulated by this contradiction, we have 
examined the conditions under which the hydrodynamic flow in the interior of a self- 
gravitating perfect fluid can become dynamically equivalent to geodesic motion in the 
same source. By definition, dynamical equivalence rests on the functional similarity 
between the corresponding differential equations of motion. In this case, the spaces 
of the solutions of these two kinds of motion are isomorphic. In other words, given 
a solution to a problem of hydrodynamic flow in a gravitating perfect fluid, one can 
always construct a solution formally equivalent to the problem of geodesic motion in 
the same source and vice versa. 

As regards a self-gravitating perfect fluid with metric g m i, dynamical equivalence 
between the two kinds of motion can be achieved in the case of adiabatic flow. Then, 
we actually transfer the problem to a virtual fluid, with metric g m i, in which the 
motions are geodesies. In this case, the two metrics are connected by means of 
a conformal transformation [Eq. (13)], with the identity element representing an 
isobaric flow of the original fluid. 

In general, to determine the functional form of the conformal factor Q(x l ) , involves 
the solution of the compatibility equation (24). In the present article we have consid- 
ered the particular solution $ = constant. The resulting conformal transformation 
affects both the geometry and the stress-energy tensor (T m i) of the original fluid, 
in a way that the functional form of the Einstein field equations remains invariant. 
In this respect, we have also determined the form of the virtual stress- energy tensor 
(T m i) in terms of the original one. In comparison to the measurements relevant to the 
observational data from the central region of the AGNs, these two quantities have a 
clear physical interpretation: 7~ m i is what we actually measure by assuming geodesic 
motions, while T m i corresponds to the "real" stress-energy tensor associated to the 
fluid matter content in those regions. In this case, with the aid of the weak energy 
condition, we have found that the observationally determined quantity (T m \) is being 



18 



overestimated with respect to T m i. 

Admiting an idealized model in which a spherically symmetric perfect fluid sur- 
rounds a massive black hole, we have applied the previous results in the determination 
of the central nuclear mass, as regards some particular AGNs. In the Newtonian limit, 
the corresponding results indicate that the mass-excess (m*) arising in the observa- 
tional determination of the central masses by ignoring the contribution of the fluid's 
internal characteristics, is not always negligible compared to the mass of the black 
hole. In fact, it can account from a few hundredths of thousanths (in NGC 4258) 
to several hundredths (in NGC 4261) of Mbh, depending on the linear dimension 
of the circumnuclear gas. On the other hand, we have shown that m 8 can be either 
comparable or even larger than the total rest-mass of the fluid involved. 

However, the assumption of a (spherically symmetric) perfect fluid can be rather 
naive in realistic situations. In this case, the viscous hydromagnetic flow in the in- 
terior of a warped accretion disk should be more appropriate, but the corresponding 
equations are too complicated for astrophysical applications. Furthermore, it can 
be verified that in the case of a magnetized self- gravitating perfect fluid, the deter- 
mination of the corresponding conformal transformation involves the solution of a 
compatibility condition of the form 

Vt m e c 

where, m e is the mass of the charge e and T™ is the antisymmetric tensor of the 
electromagnetic field involved. Clearly, Eq. (80) is far more complicated than Eq. 
(24). 

Finally, one cannot help asking whether dynamical equivalence between the two 
kinds of motion implies physical equivalence, in general. Clearly, the geodesic motion 
of a test-particle is an approximate description with certain limits of validity. The 
test-particle does not react back to modify the original gravitational field. To decide 
on how realistic is that concept, is an old and very interesting problem. In this 
respect, the conformal transformation (30) provides a direct way of linking the laws 
of motion of an idealized point-particle with those determining the hydrodynamic flow 
of a realistic volume element in the interior of a continuous source. In the latter case, 
the problem of backreaction is directly involved in the structure of the corresponding 
stress-energy tensor. To which extend is this procedure applicable to any form of T m i, 
remains an open problem. 
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